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ABSTRACT:  Several  states  of  septet  and  quintet  spin  of  FeO  are  studied  Bound 
excited  states  of  Z  ,  E”,  it,  and  z  symmetry  with  Rfi  o,  3.7  Bohrs  and  De 
^  1.3  eV  are  found.  It  is  concluded  that  the  ground  state  of  FeO  is  not  a 
septet  state.  It  is  possible  that  a  quintet  state  with  molecular  structure 
similar  to  that  of  the  states  studied  may  be  identified  with  the  c  state  of 
FeO. 
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I.  INTRORIcriOM 


i . 


In  the  course  of  a  series  of  Hart rcc-l;ock  calculations  on  FeO,  a  state 

of  V  symmetry  with  an  unusual  molecular  structure  was  obtained.  The 

Molecular  orbital  (MO)  configuration  of  this  state  is  -  besides  filled  shells 

2  ?  1  13  3  2 

for  the  atomic  cores  -  7o  8a"  Pa'  IOo'  3n'  4if  .16  .  'flic  unusual  aspect  of 
this  state  is  the  fact  that  the  molecular  orbitals  arc  almost  entirely  noil- 
bonding.  The  last  occupied  o  orbital,  IOo,  is  primarily  an  be  4s  orbital 
with  a  modest  amount  of  4p  hybridization.  The  4s-4p  hybrid  is  directed  out  of 
the  bond  (away  from  0).  Prom  examinations  of  the  Mil's,  particularly 
population  analyses^,  the  molecular  structure  is  seen  to  strongly  resemble 
that  of  an  Fc+  ion  plus  an  0  ion.  The  configuration  of  the  Pe+  ion  is 
3cl6(5ll)  4s ‘ 81pa * u ;  in  terms  of  orbitals  with  molecular  symmetry  it  is 


3da]  3dm3  (V)  3d57  (Y)  4s * 8  4pa*2 

7  X 

Tlic  configuration  of  0  is  2po"  2pm  . 

2  3 

The  choice  of  the  STO  basis  set  for  the  analytic  open  shell  SCP  *' 
calculations  is  discussed  in  Sec.  II  and  the  results  of  the  SCP  calculations 
on  the  l  state  are  presented  in  Sec.  III.  The  computed  equilibrium 
separation  is  3.73  Hohrs.  The  dissociation  energy  computed  with  respect 
to  the  restricted  Hartrcc-Fock  energies  of  the  neutral  ground  state  atoms 
is  1.3  eV. 

A  large  number  of  molecular  wave  functions  can  be  constructed  from 
the  model  of  Fe+  and  0  described  above.  This  is  done  by  distributing  the 
Feh  d  and  0  p  electrons  in  different  ways  among  the  nearly  degenerate 
dj,  dr,  and  dr  and  pa  and  pr  ?KVs  and  by  coupling  the  angular  momenta  of 
ihe  open  shells  in  different  ways.  Some  of  the  states  which  could  be 


constructed  in  this  way  were  investigated  and  four  pot  en t  i a  1  curves  almost 

7  t 

parallel  to  and  veiy  near  the  Y.  state  curve  were  found.  These?  slates 
were  investigated  by  means  of  very  limited  Cl  calculations.  The  results  of 
these  calculations  were  equivalent  to  fro  ten  orbital  c:  halations.  In  this 
context,  frozen  orbital  calculations  mean  calculations  in  which  the  SC!:  MO's 
of  the  T.  state  arc  used  to  construct  one  configuration  wave  functions  for 
other  states.  The  results  of  these  calculations  are  presented  in  See.  IV. 

Some  simple  arguments  are  given  in  Sec.  V  to  show  that  the  molecular 
correlation  energy  of  the  states  considered  is  very  close  to  the  correlation 
energies  of  the  separated  atoms.  It  is  concluded,  in  See.  VI,  that  the 
curves  obtained  arc  for  fairly  highly  excited  states  of  FcO.  It  is  also 
concluded  that  the  ground  state  of  1'eO  is  not  a  septet  state. 

II.  SLATl-ll  TYPH  BASIS  SFT 

A  reasonably  large  Slater  Type  (STO)  basis  set  based  on  carefully 

optimized  atomic  basis  sets  was  used  for  the  calculations  reported  in  this 

x 

paper.  As  is  well  known'  ,  an  STO  basis  function,  y>  is  defined  by: 

X  -  Nv^-hxpC-upy^CPa);  0) 

where  N  is  a  normalization  factor,  n  is  referred  to  as  the  principal 
quantum  number,  £  is  referred  to  as  the  basis  function  exponent,  and  the 
coordinates  0a,  and  $  arc  measured  with  respect  to  center  a.  The 
atomic  SCF  calculations  to  be  discussed  below  were  all  obtained  usine  an 

o 

2  7)  A 

open-shell  restricted  llartrec-FocK  *'  formalism.  In  this  context., 
restricted  Uartrcc-Fock  mums  that  all  the  orbitals  in  a  shell  arc  con¬ 
strained  to  have  the  same  radial  parts. 


5 


The  starting  point:  for  the  bar. is  i unctions  centered  on  Te  vans  an  8s, 

r, 

Sp,  and  ‘Id  SIO  basis  set  optimized  for  the  "  1)  state  of  l'e  by  Bagus  and 

t,  2  3 

Gilbert'.  This  basis-  set  yields  a  restricted  Hart vee-l'oc.k  *'  energy  of 

-126?.  ‘1-127  a. u.  for  the  ‘'p  state  of  lc.  (Clement  i*’  us  in,;*  a  basis  set:  of 

Us,  6p,  and  5d  STO's  obtained  -1262.4425  a.u.  and  tlic  energy  obtained  from 

.  4  7 

numerical  integration  of  the  llartreo-Tock  intogvo-differential  equations  ’ 
is  1262.444  a.u.)  This  atomic  basis  set  for  Pc  was  supplemented  with  2  4p 
and  2  4f  S  1'0's  for  the  molecular  PcO  calculations.  The  two  4p  exponents  were 
chosen  to  be  approximately  10?.  smaller  than  the  4s  exponents  in  the  atomic 
To  basis  set.  Those  exponents  were  added  to  allow  for  polarization  of  the 
4s  shell  of  be.  The  two  4f  exponents  were  chosen  to  span  roughly  the  same 
space  as  the  three  smaller  3d  STO's.  These  additional  exponents  were  not 
optimized. 


The  starting  point  for  the  basis  functions  centered  on  0  was  a  5s 

and  5p  STO  basis  set  optimized  for  0’  by  dementi*"’.  This  basis  yields  an 

SCb  energy  of  -74.78948  a.u.  for  0  (2p) .  ..'o  have  also  calculated  fin  SCF 

wave  function  neutral  0  with  the  same  basis  set  and  obtain  an  energy  of 

-74.8065  a.u.  for  0  (3p) ;  this  should  be  compared  with,  the  energy  of 

-74.S093S  a.u.  for  0  (3p)  obtained  by  Bagus  and  Gilbert^  by  optimizing  the 

4  7 

basis  set  for  neutral  0.  The  numerical  Ilartree-Pock  ’  energies  for  0  (3p) 

and  0  (2p)  arc,  respectively,  -74.80941  a.u.  and  -74.7897  a.u.  Thus  this 

0  basis  set  is  quite  satisfactory  for  calculations  on  atomic  0  and  O’. 

This  basis  set  was  supplemented  with  3  3d' STO’s  for  the  molecular  FeO 

calculations.  The  3d  exponents  are  the  same  as  have  been  used  in  several 

8 

recent  calculations  on  other  metal  monoxides  (LiO,  A£0,  MgO  and  TiO)  and 


CO,  NO,  and  0?. 


The  ban. ir.  set.  for  the  molecular  1'cO  ea leu] at  ions  was  formed  from  the 
one  center  STO's  on  Fe  and  0,  described  above ,  by  cons!  met  ing  all  possible 
basis  functions  of  o,  n ,  and  <!  symmetries.  The  details  of  the  basis  set 
are  given  in  Table  1.  Iron  previous  experience  with  diatomic  molecular  SCF 

C) 

calculations'  ,  it  is  reasonable  to  estimate  that  the  SCF  energies  obtained 
for  beO  using  this  basis  set  wil  l  be  within  .01  u.u.  of  tin  Hart  rco-Toch  limit  . 
Ill.  SCI-  Ul7, SLUTS  TOR  nn:  V  STATJ;  (\!:  TeO 


The  SCI7  calculations  reported  here  were  pcrfontied  using  the  ALG1F.MY 
quant  ion  chemistry  program  system^.  The  open  shell  SCF  fojmnlism  used  is  an 

7 

extension  of  Roothaan's  analytic  open  shell  SC!  formalism'  to  systems  with 

3 

more  than  one  open  shell  per  symmetry.  As  with  Roothaan's  earl  ter  lomnlism'  , 
the  matrix  SCF  equations,  including  the  off-diagonal  Lagrange  multipliers, 
arc  solved  without  approximation. 


The  electronic  configuration  of  the 


7„+ 


st  ate 


la2  2a2  3o2  4a2  5a2  6a2  7 a2  8a2  9a1  10a1  1H4  2II4  3H3  4 II*'  1£2. 


The  la  through  6a  and  in  and  2JT  110's  represent  the  atomic  cores; 
i.c.,  the  Is,  2s,  2p,  3s,  and  3p  orbitals  of  1'c  and  tbc  Is  orbital  of  0. 
No  further  reference  will  be  made  to  these  orbitals.  The  coupling  of  the 
remaining  shells  is  as  follows: 


I,  20  2 

i  /a  8a 


7  +  7  +  ^  + 

9a  CZ  )10c(  Z  );  L 


3  ?  3  ?  3  _ 

3ir’S(Oi)4ii'5(  II) ;  h 


,►2^3.,*^)  ,  7.,+ 
16  (  L  )j  ;  L  , 


The  SCF  wave  function  for  the-  configuration  described  above  is  invariant 
under  a  unitary  transformation  of  the  9a  and  10a  orbitals.  A  unique 
choice  of  these  orbitals  was  made  by  setting  the  off-diagonal  Lagrange 


multipliers  e^.  ()  =  0.  This  is  equivalent  to  Die  usual  choice 

’  '  ]  1 

of  cannon i cal  SCF  orbitals  for  closed  shells  . 

SCI*  calculations  linve  been  performed  at  6  internucl car  separations: 

R  =  3.2  (.2)  4.2;  the  computed  potential  cm ve  is  given  in  figure  1.  In 
Table  II,  v.’o  give  population  analyses^,  expectation  values  of  z  with 
respect  to  !’e  and  0,  and  orbital  energies  for  the  7o  tli rough  10a,  Sir, 

4  i,  and  15  orbitals  for  each  computed  internuclcar  separation,  be  also  give 
total  energies,  values  of  Y/T,  dipole  moments,  and  the  gross  atomic  popula¬ 
tions.  The  r.-nxis  is  the  internuclcar  axis  and  is  directed  from  Pc  to  0. 

The  calculated  dissociation  energy  and  equilibrium  separation  for  the 
1 state  are  =  1.2  cV  and  R ^  -  3.73  a.u.  The  dissociation  energy  is 
computed  with  respect  to  the  restricted  Hart  roe -Pock  energies  of  the  neutral 
ground  state  Pc  and  0  atoms.  These  energies^  arc  H  [0 (^P) ]  =  74.8094  a.u. 


8  .  • 


and  R[  be  (’!>)] 


1202.4427  a.u.  The  minimum  of  tb.o  SCF  potential  curve  for 


FeOwas  determined  by  making  a  quadratic  fit  to  the  calculated  energies  for 

R  =  3.0,  3.8,  and  4.0  a.u.  A  Dunlmm  analysis1^  of  the  potential  curve  gave 

gj  =  589  cm  1  and  w  x  =  3  cm  ^  for  ^0  and  '^Fe. 
e  e  e 

Although  there  are  small  variations  along  the  curve,  the  features  of  the 
ID's  arc  very  much  the  same  for  all  the  values  of  R  considered.  From  the 
population  analyses  and  expectation  values  of  z  given  in  Table  II,  we 
draw  the  following  conclusions.  The  7 a  orbital  is  almost  entirely  0  2s 
with  a  small  mixing  of  basis  functions  on  Pc.  The  8a  orbital  is  'vPCDo  0  2po 
with- alO°;  of  the  charge  on  Pc;  the  charge  on  Pc  decreases  as  the  internuclcar 
separation  decreases.  The  9a  orbital  is  almost  entirely  3da  on  Pe.  The  10a 
orbital  is  c80°d  Pc  4s  and  ^20°  Pc  4po;  the  negative  value  of  <z>  with  respect 
to  Fe  shows  that  the  sp  hybridization  is  such  that  the  charge  is  directed 


V . 


6 


away  from  the  bond.  The.  3tt,  4  r  and  16  orbitals  are  non -bonding  atomic 
orbitals  which  are  respectively  3dir  on  he,  2pir  on  0  and  3d 6  on  be. 

In  summary ,  the  charge  distribution  is  nearly  that  of  R  and  0  ; 
where  for  l’e+  tire  valence  shell  structure  (including  the  3d  .shell)  is 
3do^  3cfi;^  3d£ ^  4s  *  ^  4ps '  ^  (S)  and  for  0  ,  it  is  2p.  ^  2] jtt'^  (*41).  ’Ihe 
d- shell  electrons  of  Pc’1  arc  coupled  to  ^1)  which  is  by  ! kind's  rule  the 
energetically  most  favorable  coupling,  be  note  also  that  the  wave  function 
lias  the  flexibility  to  dissociate  to  the  ground  states  of  flic  separated 
be  and  0  atoms . 

Decks  of  cards  summarizing  the  SCI7  results,  and  including,  in 
particular,  the  expansion  co-efficients  for  the  MO's,  have  been  prepared. 
There  arc  158  cards  for  each  intemuclenr  separation.  Copies  of  these 
decks  are  available  upon  request,  from  the  author. 

IV.  FROZEN  ORBITAL  RESULTS  LOR  STATES  RELATED  TO  11  ih  V  STATE  Of  PcO. 

If  an  atomic  model  of  be+  (sit)1  and  0  p^  is  assumed,  there  are  a 
large  number  of  ways  in  which  the  Fc  3d  electrons  can  be  distributed  among 
da,  dm,  and  d5  and  the  0  2p  electrons  among  pa  and  pm.  Even  if  wc  restrict 
ourselves  to  cases  where  the  6  d-electrcns  are  coupled  to  ‘\l,  there  are 
still  a  large  number  of  septet,  quintet,  and  triplet  states  which  can  be 


constructed.  In  Table  III,  wc  list  all  possible  septet,  states  and  the 

5  + 

E  states  which  can  he  constructed  on  this  model. 


7.,+ 


Wo  have  used  the  'e'  SCF  orbitals  "  to  construct  wave  functions  for 
the  ^E+,  ^E  ,  'n,  and  ^E+  configurations  listed  in  Table  III.  For  ^E+, 
we  have  also  constructed  tlic  nine  additional  configurations  which  arise  when 
the  Fe  3d  electrons  are  not  restricted  to  be  coupled  to  ^D.  Clearly,  no 


additional  septet  configurations  can  be  constructed  by  dropping  the 
restriction  to  ^D  coupling  for  be  d^.  Configuration  interaction,  Cl 


7. 


calculations  were  then  performed  using  these  wave  functions. 

The  off-diagonal  Hamiltonian  matrix  elements  among  the  configurations 

S  + 

listed  in  Table  III  (and,  in  the  case  of  'X  ,  between  the  four  configurations 
listed  and  the  nine  addition  il  ones  mentioned  above)  are  negligible.  The 
Cl  eigenfunctions  could  alwavs  bo  clcarlv  identified  with  a  single  one  of 
the  configurations  listed  in  Table  111;  the  dominant  Cl  coefficient  for 
those  stales  is  always  greater  than  0.905.  Thus  the  results  for  these  states 
are  the  same  as  would  be  obtained  by  frozen  orbital  calculations.  (This 
was  not  expected  before  the  CT  calculations  were  performed  but  follows  for 
the  reasons  given  next.)  The  off-diagonal  matrix  elements  arc  small  because 
of  the  atomic  nature  of  the  Mil's  and  the  sorts  of  replacements  allowed  between 


the  configurations.  In  brief,  these  off-diagonal  matrix  elements  involve 
only  two  center  integrals  between  MO's  which  are  essentially  centered  on  be 
or  on  0.  The  exchange -type  integrals,  [FcO|FcO],  arc  small  because  the  over¬ 
lap  between  the  Fc  and  0  orbitals  is  small.  The  Coulomb-type  integrals, 
[lcFe|0  0],  are  small  because  the  angular  factors  involved  are  such  that  the 

integrals  represent  high  order  multipole  interactions.  There  are  large 

5  + 

interactions  among  some  of  the  additional  L  configurations.  However,  the 

lowest  CT  eigenvector  significantly  involving  one  of  these  configurations 

c  +  5  + 

is  the  forth  not  of  the  T,  CT  and  is  ^2.5  cV  above  the  lowest  l  root; 

thus  these  states  arc  not  of  particular  interest. 

7  +  7  -  7 

The  calculated  potential  curves  for  all  the  T,  ,  l  ,  and  II  states 
and  for  the  lowest  two  of  the  states  are  given  in  Figure  1.  'Die  Cl 
results  for  the  lower  state  (configuration  1  of  Table  III)  are  the  same 


as  the  SCF  results  to  eight  significant  figures.  The  curves  for  the  remain- 
5  + 

mg  eleven  ‘7  states  are  not  bound  with  respect  to  separated  Hartree-Fock 
atoms  and  are  above  the  scale  of  Figure  1.  In  Table  IYr,  we  tabulate  the 


p 


calculated  R  and  I)  for  each  of  the  states  plotted  in  Figure  J.  The 
c  c 

com  inant  coni  i  gurat j  on  in  each  of  the  .states  is  also  identified  by  refe  rcnce 

to  Table  ITT.  The  values  of  R  and  1)(  liavc  been  oiitaivieci  by  fit. tiny,  a 

quadratic  around  the  minimum  of  the-  curves.  As  usual,  b  is  defined  with 

respect  to  the  llartrce-Fook  energies  of  the  separated  atoms. 

For  cadi  symmetry ,  the  lowest  state  lias  1^  'u  3.7  a.u.  and  D  ^  1.2  - 

1.5  eV.  Tt  is  quite  likely  that  there  will  be  states  with  values  of  R 

and  I)  close  to  those  found  here  for  the  other  septet  and  quintet  symmetries 

7  7  5  -  5  5  5 

(  A,  ]),  'Z  ,  ‘n,  ‘A,  and  *  <ji)  and  possibly  also  for  some  of  the  triplet,  states 

arising  from  Fe*  [ cl ^ ("'’!))  (sp)]  and  0  . 

Of  course,  a  variational  treatment  of  the  states  listed  in  Table  IV 
is  likely  to  change  their  energies.  However,  the  point  of  these  frozen 
orbital  calculations  is  to  show  that,  there  arc  several  states  of  different 
symmetries  but  similar  R  end  D  .  Because  the  binding  is  almost  entirely 
ionic  and  the  ID’s  essentially  atomic  orbitals,  it  is  not  likely  that  the 
frozen  orbital  results  would  be  changed  greatly  by  variational  calculations. 
V.  ESTIMATES  OP  THE  iDLECULAR  EXTRA  CORRELATION  ENERGY. 


The  Hartree-Fock  and  experimental  transition  energies  and  the 
correlation  energy  changes  for  transitions  from  atomic  Fe  to  PeH  and  0  to 
0  are  tabulated  in  Table  V.  The  correlation  energy  of  Fc  (  D)  +  0  (“p) 
is  increased  by  +0,38  cV  with  respect  to  the  neutral  atoms.  If  we  follow  the 
population  analysis  given  in  Table  II  and  assume  that  Fc+  in  FeO  is  vSOp 
Fe+  3d  Ms (^D)  and  v20o  Fc+  3d(l4p(",l:)  then' the  correlation  energy  is  increased 
by  only  +0.23  eV.  On  the  model  that  the  states  of  FeO  considered  here  are 
basically  separated  ions;  t hi s  means  that  the  amount  of  extra  molecular 
correlation  energy  is  small  and  about:  .2  -  .4  cV.  Thus  the  computed 


1 


dissoci ;>i  Lou  energies  given  in  Table  IV  are  libel)'  to  be  reasonably  close 
to  the  true  values. 

V.  CONCLUSIONS 

Harrow  and  Senior^’  have  suggested  that  the  mast  likely  ground  state 
7  r) 

for  l?cO  is  a  Y,  or  '  Y,,  We  believe  that  we  have  investigated  tlie  possible 

septet  configurations  fairly  carefully  and  it  seems  \ery  unlikely  that  these 

dissociation  energies  can  be  significantly  larger  than  those  given  in  Table  IV. 

brewer  and  Rosenblatt1 '  have  estimated  a  value  of  D  °  -  9.1  ±  kcal/mole 

(4.12  ±  .22  cV)  for  ReO.  A  veiy  recent  t lie roio- chemical  value  of 

Do°  =  97  ±  3  kcal/mole  (4.21  i  .13  cV)  has  been  reported  by  Rnlducei.  et  al w. 

7 

Tills  seems  to  rule  out  a  septet  state  and  in  particular  a  T.  state  as  the 
ground  state  of  TcO. 

Several  quintet  states  with  structure  quite  different  from  the  states 

considered  here  have  been  investigated  and  a  large  Cl  calculation  performed 

on  a  ‘V,-1  state.  lor  the  present,  we  note  that  assuming  a  dissociation  energy 

of  'v  4.2  eV,  tlien  tlie  states  considered  in  this  paper  will  have  term  energies 

19 

of  'v  2.7  eV.  'Ilierc  are  observed  states  of  Fed  ‘  ,  the  a,  b,  and  c  states 


with  term  energies  cf,  respectively,  2.14,  2.22,  and  2.79  cV.  If  the  ground 

state  of  FeO  is  a  quintet  state,  then  it  is  possible  to  identify  the  c  state 

with  a  quintet  state  with  the  structure  which  we  have  considered  in  this 

19  -1 

paper.  Tlie  observed  '  w  -  S40  cm  for  tlie  c  state  is  in  keeping  with 
the  shallow  curves  observed  for  the  states  reported  here. 
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TAI1LL  IV. 

Values  of  Re  and  He  for  Frozen  Orbital  Calculations 
on  FcO.  The  dominant  conf iquration  for  encli  state 
is  identified  by  reference  to  the  list  qiven  in 
Table  111. 


Symmetry 

Dominant 

Configuration 

Re(a.u.) 

n0(ev) 

V  * 

1 

3.73 

1  .24 

V 

2 

3.90 

-0.60 

V 

3 

3.74 

1  .21 

h 

5 

3.G9 

1  .46 

7n 

6 

3.83 

0.79 

7« 

4 

3.72 

0.06 

V' 

10 

3.74 

1.15 

V 

11 

3.71 

0.20 

*Those  results  arc 
ci ah t  sionificant 

identical  to  the  SCF 
f i oures . 

results  for  this 

state  to 

18. 


TAB I T  V. 

Calculated  and  rx,rtf;r'i r.i^r.tc  1  Trorr  i  lion 
Enei  cii os  and  Correlation  I  ncrgy  Chaiv'  s 
in  Fe  and  0.  Enemies  arc-  in  eV. 


State 

Calculated 

HF  Energies 

Exneriii.oiv!  al 
Enon-  i :  s 

AE  Correlation 

re  3d6 (^D) 4s 2 

0 

0 

- — 

Fe+3df,(bD)4s;  S 

6.275 

7,891)  b 

-1.623 

Fe+3dVo)4p;  6F 

10.751 

13.101  h 

-2. 350 

0  2p4(3p) 

0 

0 

0"  2p5(2p) 

+0.541 

-1.456  c 

+2 . 006 

a. 

The  HF  energies 
using  a  program 

are  obtained  from  numerical  Hartree-Fock  calculations 
of  C.  Froese-Fischer.  See  Ref.  7. 

b. 

See  Ref.  14. 

c. 

Sec  Ref.  15. 

Energy  (A.U.) 


Figure  1.  Potential  curves  for  several  Septet  and  Quintet  States  of  FeO. 
The  lower  k  curve  is  obtained  from  SCF  calculations;  the  remaining 
curves  are  obtained  from  frozen  orbital  calculations  using  the  /'£+SCF 
WO’s  as  describ'd  in  Sec. III.  The  dissociation  limit  is  for  the  SCF  energies 
of  the  ground  state  Fc  and  0  atoms. 


